Statistics of microscopic yielding in sheared aqueous foams 



Yuhong Wang, Kapilanjan Krishan, and Michael Dennin 

Department of Physics and Astronomy, University of California at Irvine, Irvine, California 92697-4575 

(Dated: February 6, 2008) 

We detail the statistical distribution of bubble rearrangements in a sheared two-dimensional foam. 
Such rearrangements, known as Tl events, are vital to mechanisms resulting in flow through micro- 
scopic mechanical yielding. We find that at a constant rate of shear, the rate of occurrence of Tl 
events shows only small fluctuations. This rate is however seen to vary significantly with a variation 
in the initial configuration of bubbles constituting the foam. In addition, we detail the spatial and 
orientational distribution of Tl events and relate them to the distribution of stresses in the bulk of 
the material. Some insights into the irreversibility of the dynamics are also discussed. 
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I. INTRODUCTION 

In many materials, the microscopic structure predom- 
inantly determines the mechanism of flow. The micro- 
scopic response to deformation may vary from linear 
in the case of elastic materials and newtonian fluids to 
highly nonlinear in non-newtonian fluids and granular 
materials. In many materials, the response is linear for 
sufficiently small deformations, only to become nonlin- 
ear at a critical value of the deformation. The transition 
to nonlinear behavior is often controlled by the nature 
of various microscopic events. It is worth highlighting 
just a few examples. In crystalline solids, defects in the 
lattice structure may lead to fractures within the ma- 
terial l]. In amorphous materials, plasticity has been 
described in terms of local rearrangements ofparticles, 
known as shear transformation zones (STZs) For 
soft materials, such as colloids, foams, and emulsions, a 
soft glassy rheology (SGR) model based based on par- 
ticle rearrangements is able to capture many features of 
the macroscopic behavior 0, Q . Similarly, in the con- 
text of emulsions, a model based on "weak" zones within 
the material has been successful in explaining anomalous 
measurements of rheological properties ja,l3- Therefore, 
characterizing microscopic rearrangements is imperative 
to understanding the behavior of materials subjected to 
external loads during stress, strain etc. 

In this article we study the behaviour of aqueous foams 
when subjected to macroscopic shear. The yielding of 
such foams is typically described by a Bingham model 
that proposes a critical yield stress at which the mate- 
rial switches from an elastic to a plastic response [||. 
When stresses in the material exceed the yield stress, 
failure in the microscopic structure results in flow within 
the material. The microscopic aspect related to this be- 
haviour is a succession of slip between the constituent 
bubbles of the foam, known as Tl events (see Fig. 
In its most basic manifestation, such slip involves four 
neighboring bubbles, with one pair of nearest neighbours 
becoming next-nearest neighbours and vice- versa for an- 
other pair. Tl events have been the subject of study 
in simulations 0, 0, 0> 0- indirect measurement in 
three-dimensional foam using diffusive wave spectroscopy 




FIG. 1: The above illustrates the rearrangement dynamics 
during a Tl event. The centers of complimentary bubble 
pairs that undergo neighbor switching are joined by line seg- 
ments. The two vertically connected bubbles become neigh- 
bours^) from being next-nearest neighbours(NN), symbol- 
ized by NN— >N, and vice-versa for the horizontal bubbles, 
denoted by NN^N. 

In addition to shear induced Tl events, flowing foams 
are subject to a number of other phenomenon such as 
drainage, aging, coarsening and T2 events (the disap- 
pearance of a bubble within a foam) 8]. While these 
mechanisms do have an influence on the flow behaviour, 
their role is minimal at the short time-scales of response 
we are interested in. Various other studies look at the 
interaction of these mechanisms with flow, for example 
drainage- induced convection |21| . 

The role of Tl events resulting in flow in aqueous foams 
has a number of interesting implications. The spatial 
distribution of these events would provide clues to the 
distribution of internal stresses, as Tl events would be 
expected to occur more frequently in regions where the 
stresses exceed the local yield stress of the material. A 
similar argument may be used in regard to the orienta- 
tion of the Tl events. The slippage between neighbouring 
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bubbles along a preferred direction indicates an orienta- 
tional inhomogeneity within the material or the imposed 
forces. The temporal distribution of Tl rates would also 
be insightful into relaxation processes that occur: A burst 
of Tl events between periods of relative calm could indi- 
cate avalanche- type behaviour wherein stresses are built 
up slowly over time and released rapidly 0, l22| . And fi- 
nally, a relationship between the configuration of bubbles 
and the spatial and temporal distribution of Tl events 
help build connections between the structural and dy- 
namical properties of such materials. 

In investigating these areas, it is important to develop 
a large statistics detailing the spatial and temporal dis- 
tribution of Tl events. Though DWS d Q Hi E E3 
has been very successful in providing statistical infor- 
mation on individual bubble motions, there are no well- 
established techniques to gather data within the bulk of 
a three-dimensional foam at the time-scales and spatial 
resolution we require to track Tl events during shear. We 
therefore focus our attention on a two-dimensional foam, 
also known as a bubble raft. The bubble raft has been 
found to be a useful guide in the study of foams 0, |2^| 
as well as solid materials 0, 0, H3- Also, the impact 
of boundary conditions on Tl events in bubble rafts has 
been characterized 20]. In addition, it is helpful to ini- 
tially study monodisperse bubble rafts. This makes it 
easier to interpret the results and eliminates complexi- 
ties associated with polydispersity . Such effects may 
include faster coarsening rates, the development of spa- 
tial inhomogenities due to size dependent flow behaviour 
of individual bubbles etc. We also focus our observations 
in a regime when the foam is subjected to high rates of 
shear. It is expected that the system is not in a qua- 
sistatic regime, and the response is far from equilibrium. 
In contrast, prior studies have detailed the response at 
lower rates of strain [j| 0, ^ , and the creep response 
of foam in which the impact of coarsening induced Tl 
events is studied [l^. 

In the next section we detail techniques of our ex- 
periments as well as image analysis. Following this, we 
present and discuss our results on the distribution of Tl 
events under steady shear. 



II. EXPERIMENT AND IMAGE ANALYSIS 

Our experiments utilize a layer of monodisperse bub- 
bles floating on the surface of water undergoing shear 
between two counter- rotating bands. The bubbles were 
made in a trough containing an aqueous solution con- 
sisting of water, glycerol and miracle bubble solution in 
the ratio of 80%, 15% and 5% by volume respectively. A 
small amount of fluorescent dye, Fluorescein, was added 
as well to enhance imaging within our system. A steady 
stream of nitrogen injected into this solution using a nee- 
dle was sufficient to generate a monodisperse layer of bub- 
bles, resulting in a bubble raft. 

In addition, two bands with evenly spaced gaps along 



their length were positioned at the surface of the bubble 
monolayer with gears attached to submerged shafts. The 
shafts are supported by teflon bearings attached to the 
bottom of the trough. The bands were attached to a step- 
per motor through which their rate of motion could be 
accurately controlled. The bubbles floating in between 
these bands undergo linear shear. The amount of shear 
is quantified by 7 = 2L(t)/D, where D(= 6cm) is the 
width between the bands and L the distance traversed 
by the band in time t. Typically, the width of the band 
corresponds to about 20 bubble diameters. The entire 
trough was enclosed, and illuminated using UV lamps. 
The motion of the bubble raft was recorded as consecu- 
tively saved images using a CCD array. The frame rate 
of image capture was maintained high enough so as to 
identify individual bubbles between successive images. 

The images captured were filtered to remove noise in- 
herent in the CCD imaging technique as well as some 
optical inhomogenities associated with the optics. The 
filtered images were thresholded to yield a binary image 
delineating the location of individual bubbles. The center 
positions of individual bubbles were computed as the cen- 
ter of mass of each of the delineated regions. A voronoi- 
construction was then used to construct the edges be- 
tween the bubbles constituting the bubble raft. Nearest- 
neighbours were designated as those that shared a com- 
mon edge. Next-nearest neighbors were also extracted 
using this information. Individual bubbles were identi- 
fied based on the least displacement of bubbles between 
consecutive images. 

The identification of Nearest-neighbors and next- 
nearest-neighbors between successive frames was used to 
establish the occurrence of Tl events: i.e. when a set of 
of nearest neighbors became next nearest neighbors. The 
location of the Tl event was designated as the center of 
mass of the bubbles participating in the event. A more 
detailed account of the experimental and image analysis 
techniques may be found in reference |26J. The exper- 
iments are performed over a number of rates of shear, 
7 = dry/ dt, producing a sequence of images taken at 
equally spaced time intervals. 



III. RESULTS AND DISCUSSION 
A. Spatial distribution: Slip zones, homogeneity 

The Tl events are seen to be distributed uniformly 
between the bands in the region of interest. There is a 
slight increase in the number of Tl events in the vicinity 
of the bands. This may be attributed to pinning effects at 
the bands resulting in locally higher shear stresses. The 
spatially uniform distribution of Tl events seen in the 
lower plot in Fig. [2] suggests that the material properties 
as well as forces imposed are homogeneous when averaged 
over the entire period of shear. 

Often, multiple Tl events occurs along a local crys- 
talline axis of the foam as in the upper plot in Fig. [21 
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FIG. 2: The above two plots indicate the positions of Tl 
events over different strains. The lower figure is strained 
about ten times longer than the upper figure. The top fig- 
ure indicates the succession of Tl events that lead to slip. 
The homogeneous spread of Tl events is seen in the lower 
figure. The increase in Tl events close to the driving bands 
(towards the top and bottom) are attributed to the pinning 
of bubbles at these boundaries. 





























FIG. 3: The blue dots in the above plot indicate the instan- 
taneous locations of the centers of bubbles forming the bub- 
ble raft. The red lines connect bubbles that have undergone 
compression, becoming nearest neighbors from being next 
nearest neighbors (NN— »N). The green lines connect bubbles 
that have undergone dilation, becoming next-nearest neigh- 
bors from being nearest neighbors initially (NN<— N). The ori- 
entation of these two events is measured by the angle (6) the 
red and blue lines make with the x-axis. The horizontal ar- 
rows at the top and bottom of figure show the direction of 
shear imposed by the bands. 



The spatial extent of such regions may extend across the 
entire width of the sheared foam, indicating the role of 
Tl events in causing macroscopic slip and flow. Such 
coherent rearrangements over small 7 implies that the 
local crystalline structure may influence the focussing of 
stresses as well as provide principal directions for fracture 
over regions of the order of the system size. 



B. Orientational distribution: Anisotropy, 
irreversibility 

In order to study the orientation of Tl events, we 
look at a single pair of bubbles that participate in the 
Tl event. This pair represents the bubbles that be- 
come nearest neighbours (N) from being next-nearest- 
neighbours (NN). The geometry used to define the an- 
gles of the Tl events is illustrated in Fig. |3| We focus 
on the angle between the x-axis (orientation of the driv- 
ing bands) and the line joining nearest-neighbors as one 
parametrization of the orientation of the Tl events. This 
is plotted in Fig. 0Jt. We find that the orientation of Tl 
events shows an distribution peaked about a particular 
angle, which may indicate a principal shear axis. At lower 
rates of shear, a secondary peak is also seen. This sug- 
gests the existence of a secondary axis for shear during 
Tl events. 

A similar measurement on the complimentary pairs of 
bubbles is shown in Fig. 0Jd. The angles corresponding 
to the maxima in the curves indicate the most preferred 
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FIG. 4: The upper figure plots the angle with respect to the 
driving bands at which next-nearest bubbles become near- 
est bubbles(NN— >N). The lower figure indicates the compli- 
mentary mechanism of NN<— N. The different colours indicate 
plots obtained at different rates of shear. 



direction of yielding through Tl events through compres- 
sion and dilation. Distinct orientations for bubbles com- 
ing closer (NN— >N) and bubbles moving apart (NN^N) 



4 



suggest an anisotropy in the internal stresses on the bub- 
bles. The difference between the maxima of angles seen 
in the upper and lower plots of Fig. 0] correspond to the 
angle between the lines indicated in Fig. ^ further il- 
lustrating the role of local structure in directing internal 
stresses. 

In addition to the difference in most probable angular 
orientation, there is also a noticeable difference in the 
widths of the distributions of Fig. The upper plot 
(NN— >N) is observed to be broader and lower than the 
lower (NN<— N). This has important ramifications on the 
reversibility of Tl events. In the ideal case of reversibility, 
all bubbles undergoing NN^N would undergo NN^N 
on reversing the shear direction/time. In such a sce- 
nario, the width of angular distributions would remain 
unchanged between NN^N and NN^N as the two pro- 
cesses are symmetric under time reversal. Our observa- 
tion require that such a reversal would need to occur at 
different ranges of angular distributions. This makes the 
reversibility of Tl events nontrivial in a statistical sense 
while placing no constraints on the reversibility of indi- 
vidual Tl events. 



C. Temporal distribution: Rate invariance, 
extensivity 
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FIG. 5: (a)Number of Tl events vs. time. At the lowest rate 
of shear, we also include the data for longer times, illustrating 
the rate invariance over longer times. (b)Rate of Tl events 
vs. rate of shear. The solid line is a linear fit passing through 
the origin. 



We also look at the rate at which Tl events occur. We 
find that even over large strains, the rate of Tl events is 
uniform, i.e. the rate of Tl events does not vary much 
during the run. It is particularly noteworthy that the rate 
of Tl events shows only a small deviation during shear 
as indicated in Fig. [5J The consistent and well defined 
slope in upper figure has a number of implications. 

Firstly, it indicates that the rate of Tl events may be 
considered equivalent to a short time constraint on the 
dynamical evolution of the system. This implies that the 
response to strain at the boundaries imposes an internal 
consistent rate of slippage between bubbles. The rate of 
Tl events is seen to increase linearly with the imposed 
rate of shear at the boundaries. 

Second, when coupled with the homogeneous spatial 
distribution of Tl events, it is relevant to talk about the 
rate of Tl events per bubble. This implies that the rate 
of Tl events is an extensive quantity. 

IV. CONCLUSIONS AND OUTLOOK 

We have undertaken a systematic study of Tl events in 
two-dimensional sheared aqueous foams. Tl events con- 
stitute the fundamental dynamics associated with foam 
flows. The study has emphasized the role of local orien- 
tation of bubbles and the resulting anisotropy in orienta- 
tional distribution of Tl events. Future work will focus 
on connecting the observed dynamics of the Tl events 
with the concept of "weak zones" developed in the con- 
text of emulsions 0, QJ and the analysis of STZs in the 
context of model foams 0. Also, we have demonstrated 
a clear correlation between the rate of Tl events and the 
rate of strain. 

We also open many interesting avenues of investiga- 
tion. For example, contrasting the statistics obtained 
with polydisperse bubble distributions would be useful 
in understanding amorphous materials. Also, the signif- 
icance of structure at the ends of Tf slip zones in Fig. [2 
could unravel mechanisms for fracture nucleation and ter- 
mination. It would be particularly interesting to under- 
stand the origins of uneven distribution widths of angu- 
lar distribution between NN^N and NN<— N. Designing 
a shear geometry that maintains symmetry in the two 
distributions would illustrate the influence of boundaries 
in determining reversibility. Finally, developing a bet- 
ter understanding of the connection with the results for 
coarsening generated Tl events in a three-dimensional 
foam [l^ wn l be an important step. 
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